High Throughput SARS-CoV-2 Genome Sequencing from 384 Respiratory Samples Using the Illumina COVIDSeq Protocol.

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic has fostered the use of high-throughput techniques to sequence the entire severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome and track its evolution. The present study proposes a rapid and relatively less expensive sequencing protocol for 384 samples by adapting the use of an Illumina NovaSeq library to an Illumina MiSeq flow cell instrument. The SARS-CoV-2 genome sequences obtained with Illumina NovaSeq and those obtained using MiSeq instruments were compared with the objective to validate the new, modified protocol. A total of 356 (94.6%) samples yielded interpretable sequences using the modified Illumina COVIDSeq protocol, with an average coverage of 91.6%. By comparison, 357 (94.9%) samples yielded interpretable sequences with the standard COVIDSeq protocol, with an average coverage of 95.6%. Our modified COVIDSeq protocol could save 14,155 euros per run and yield results from 384 samples in 53.5 h, compared to four times 55.5 h with the standard Illumina MiSeq protocol. The modified COVIDSeq protocol thus provides high quality results comparable to those obtained with the standard COVIDSeq protocol, four times faster, while saving money.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule.

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

First cases of infection with the 21L/BA.2 Omicron variant in Marseille, France.

The SARS-CoV-2 21K/BA.1, 21L/BA.2, and BA.3 Omicron variants have recently emerged worldwide. To date, the 21L/BA.2 Omicron variant has remained very minority globally but became predominant in Denmark instead of the 21K/BA.1 variant. Here, we describe the first cases diagnosed with this variant in south-eastern France. We identified 13 cases using variant-specific qPCR and next-generation sequencing between 28/11/2021 and 31/01/2022, the first two cases being diagnosed in travelers returning from Tanzania. Overall, viral genomes displayed a mean (±standard deviation) number of 65.9 ± 2.5 (range, 61-69) nucleotide substitutions and 31.0 ± 8.3 (27-50) nucleotide deletions, resulting in 49.6 ± 2.2 (45-52) amino acid substitutions (including 28 in the spike protein) and 12.4 ± 1.1 (12-15) amino acid deletions. Phylogeny showed the distribution in three different clusters of these genomes, which were most closely related to genomes from England and South Africa, from Singapore and Nepal, or from France and Denmark. Structural predictions highlighted a significant enlargement and flattening of the surface of the 21L/BA.2 N-terminal domain of the spike protein compared to that of the 21K/BA.1 Omicron variant, which may facilitate initial viral interactions with lipid rafts. Close surveillance is needed at global, country, and center scales to monitor the incidence and clinical outcome of the 21L/BA.2 Omicron variant.

SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding.

At the time of a new and unprecedented viral pandemic, many questions are being asked about the genomic evolution of SARS-CoV-2 and the emergence of different variants, leading to therapeutic and immune evasion and survival of this genetically highly labile RNA virus. The nasopharyngeal persistence of infectious virus beyond 17 days proves its constant interaction with the human immune system and increases the intra-individual mutational possibilities. We performed a prospective high-throughput sequencing study (ARTIC Nanopore) of SARS-CoV-2 from so-called "persistent" patients, comparing them with a non-persistent population, and analyzing the quasi-species present in a single sample at time t. Global intra-individual variability in persistent patients was found to be higher than in controls (mean 5.3%, Standard deviation 0.9 versus 4.6% SD 0.3, respectively, p < 0.001). In the detailed analysis, we found a greater difference between persistent and non-persistent patients with non-severe COVID 19, and between the two groups infected with clade 20A. Furthermore, we found minority N501Y and P681H mutation clouds in all patients, with no significant differences found both groups. The question of the SARS-CoV-2 viral variants' genesis remains to be further investigated, with the need to prevent new viral propagations and their consequences, and quasi-species analysis could be an important key to watch out.

Simultaneous SARS-CoV-2 Genome Sequencing of 384 Samples on an Illumina MiSeq Instrument through Protocol Optimization.

In the present study, we propose a high-throughput sequencing protocol using aNextera XT Library DNA kit on an Illumina MiSeq instrument. We made major modifications to this library preparation in order to multiplex 384 samples in a single Illumina flow cell. To validate our protocol, we compared the sequences obtained with the modified Illumina protocol to those obtained with the GridION Nanopore protocol. For the modified Illumina protocol, our results showed that 94.9% (357/376) of the sequences were interpretable, with a viral genome coverage between 50.5% and 99.9% and an average depth of 421×. For the GridION Nanopore protocol, 94.6% (356/376) of the sequences were interpretable, with a viral genome coverage between 7.0% and 98.6% and an average depth of 2123×. The modified Illumina protocol allows for gaining EUR 4744 and returning results of 384 samples in 53.5 h versus four times 55.5 h with the standard Illumina protocol. Our modified MiSeq protocol yields similar genome sequence data as the GridION Nanopore protocol and has the advantage of being able to handle four times more samples simultaneously and hence is much less expensive.

A 21L/BA.2-21K/BA.1 "MixOmicron" SARS-CoV-2 hybrid undetected by qPCR that screen for variant in routine diagnosis.

Among the multiple SARS-CoV-2 variants identified since summer 2020, several have co-circulated, creating opportunities for coinfections and potentially genetic recombinations that are common in coronaviruses. Viral recombinants are indeed beginning to be reported more frequently. Here, we describe a new SARS-CoV-2 recombinant genome that is mostly that of a Omicron 21L/BA.2 variant but with a 3' tip originating from a Omicron 21K/BA.1 variant. Two such genomes were obtained in our institute from adults sampled in February 2022 in university hospitals of Marseille, southern France, by next-generation sequencing carried out with the Illumina or Nanopore technologies. The recombination site was located between nucleotides 26,858-27,382. In the two genomic assemblies, mean sequencing depth at mutation-harboring positions was 271 and 1362 reads and mean prevalence of the majoritary nucleotide was 99.3 ± 2.2% and 98.8 ± 1.6%, respectively. Phylogeny generated trees with slightly different topologies according to whether genomes analyzed were depleted or not of the 3' tip. This 3' terminal end brought in the Omicron 21L/BA.2 genome a short transposable element of 41 nucleotides named S2m that is present in most SARS-CoV-2 except a few variants among which the Omicron 21L/BA.2 variant and may be involved in virulence. Importantly, this recombinant is not detected by currently used qPCR that screen for variants in routine diagnosis. The present observation emphasizes the need to survey closely the genetic pathways of SARS-CoV-2 variability by whole genome sequencing, and it could contribute to gain a better understanding of factors that lead to observed differences between epidemic potentials of the different variants.

Screening and Whole Genome Sequencing of SARS-CoV-2 Circulating During the First Three Waves of the COVID-19 Pandemic in Libreville and the Haut-Ogooué Province in Gabon.

Since the onset of the COVID-19 pandemic, the SARS-CoV-2 viral dynamics in Africa have been less documented than on other continents. In Gabon, a Central African country, a total number of 37,511 cases of COVID-19 and 281 deaths have been reported as of December 8, 2021. After the first COVID-19 case was reported on March 12, 2020, in the capital Libreville, the country experienced two successive waves. The first one, occurred in March 2020 to August 2020, and the second one in January 2021 to May 2021. The third wave began in September 2021 and ended in November 2021. In order to reduce the data gap regarding the dynamics of SARS-CoV-2 in Central Africa, we performed a retrospective genotyping study using 1,006 samples collected from COVID-19 patients in Gabon from 2020 to 2021. Using SARS-CoV-2 variant screening by Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) and whole genome sequencing (WGS), we genotyped 809 SARS-CoV-2 samples through qRT-PCR and identified to generated 291 new genomes. It allowed us to describe specific mutations and changes in the SARS-CoV-2 variants in Gabon. The qRT-PCR screening of 809 positive samples from March 2020 to September 2021 showed that 119 SARS-CoV-2 samples (14.7%) were classified as VOC Alpha (Pangolin lineage B.1.1.7), one (0.1%) was a VOC Beta (B.1.351), and 198 (24.5 %) were VOC Delta (B.1.617.2), while 491 samples (60.7%) remained negative for the variants sought. The B1.1 variant was predominant during the first wave while the VOC Alpha dominated the second wave. The B1.617.2 Delta variant is currently the dominant variant of the third wave. Similarly, the analysis of the 291 genome sequences indicated that the dominant variant during the first wave was lineage B.1.1, while the dominant variants of the second wave were lineages B.1.1.7 (50.6%) and B.1.1.318 (36.4%). The third wave started with the circulation of the Delta variant (B.1.617). Finally, we compared these results to the SARS-CoV-2 sequences reported in other African, European, American and Asian countries. Sequences of Gabonese SARS-CoV-2 strains presented the highest similarities with those of France, Belgium and neighboring countries of Central Africa, as well as West Africa.

Outcomes of 2111 COVID-19 Hospitalized Patients Treated with Hydroxychloroquine/Azithromycin and Other Regimens in Marseille, France, 2020: A Monocentric Retrospective Analysis.

Objectives: We evaluated the 6-week mortality of SARS-CoV-2 hospitalized patients treated using a standardized protocol in 2020 in Marseille, France. Methods: A retrospective monocentric cohort study was conducted in the standard hospital wards at the Institut Hospitalo-Universitaire Méditerranée Infection, between March and December 2020 in adults with SARS-CoV-2 PCR-proven infection. Results: Of the 2111 hospitalized patients (median age, 67 [IQR 55-79] years; 1154 [54.7%] men), 271 were transferred to the intensive care unit (12.8%) and 239 died (11.3%; the mean age of patients who died was 81.2 (±9.9)). Treatment with hydroxychloroquine plus azithromycin (HCQ-AZ), used in 1270 patients, was an independent protective factor against death (0.68 [0.52 - 0.88]). This effect was consistent for all subgroups of age, comorbidities, severity of the disease and comedications with zinc or corticosteroids. Zinc was independently protective against death (0.39 [0.23 - 0.67]), in a subgroup analysis of patients treated with HCQ-AZ without dexamethasone. The use of high-flow oxygen therapy in elderly patients who were not eligible for intensive care unit transfer saved 19 patients (33.9%). Conclusions: In our 2020 cohort, treating COVID-19 with HCQ-AZ was associated with lower mortality. These results need to be analyzed in the context of academic discussions about observational studies versus randomized clinical trials. More data will deserve to be analyzed in the SARS-Cov 2 variants, vaccination and post-vaccination era.

Sequential Appearance and Isolation of a SARS-CoV-2 Recombinant between Two Major SARS-CoV-2 Variants in a Chronically Infected Immunocompromised Patient.

Genetic recombination is a major evolutionary mechanism among RNA viruses, and it is common in coronaviruses, including those infecting humans. A few SARS-CoV-2 recombinants have been reported to date whose genome harbored combinations of mutations from different mutants or variants, but only a single patient's sample was analyzed, and the virus was not isolated. Here, we report the gradual emergence of a hybrid genome of B.1.160 and Alpha variants in a lymphoma patient chronically infected for 14 months, and we isolated the recombinant virus. The hybrid genome was obtained by next-generation sequencing, and the recombination sites were confirmed by PCR. This consisted of a parental B.1.160 backbone interspersed with two fragments, including the spike gene, from an Alpha variant. An analysis of seven sequential samples from the patient decoded the recombination steps, including the initial infection with a B.1.160 variant, then a concurrent infection with this variant and an Alpha variant, the generation of hybrid genomes, and eventually the emergence of a predominant recombinant virus isolated at the end of the patient's follow-up. This case exemplifies the recombination process of SARS-CoV-2 in real life, and it calls for intensifying the genomic surveillance in patients coinfected with different SARS-CoV-2 variants, and more generally with several RNA viruses, as this may lead to the appearance of new viruses.

Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics.

After the end of the first epidemic episode of SARS-CoV-2 infections, as cases began to rise again during the summer of 2020, we at IHU Méditerranée Infection in Marseille, France, intensified the genomic surveillance of SARS-CoV-2, and described the first viral variants. In this study, we compared the incidence curves of SARS-CoV-2-associated deaths in different countries and reported the classification of SARS-CoV-2 variants detected in our institute, as well as the kinetics and sources of the infections. We used mortality collected from a COVID-19 data repository for 221 countries. Viral variants were defined based on ≥5 hallmark mutations along the whole genome shared by ≥30 genomes. SARS-CoV-2 genotype was determined for 24,181 patients using next-generation genome and gene sequencing (in 47 and 11% of cases, respectively) or variant-specific qPCR (in 42% of cases). Sixteen variants were identified by analyzing viral genomes from 9,788 SARS-CoV-2-diagnosed patients. Our data show that since the first SARS-CoV-2 epidemic episode in Marseille, importation through travel from abroad was documented for seven of the new variants. In addition, for the B.1.160 variant of Pangolin classification (a.k.a. Marseille-4), we suspect transmission from farm minks. In conclusion, we observed that the successive epidemic peaks of SARS-CoV-2 infections are not linked to rebounds of viral genotypes that are already present but to newly introduced variants. We thus suggest that border control is the best mean of combating this type of introduction, and that intensive control of mink farms is also necessary to prevent the emergence of new variants generated in this animal reservoir.

Culture and identification of a "Deltamicron" SARS-CoV-2 in a three cases cluster in southern France.

Multiple SARS-CoV-2 variants have successively, or concomitantly spread worldwide since the summer of 2020. A few co-infections with different variants were reported and genetic recombinations, common among coronaviruses, were reported or suspected based on co-detection of signature mutations of different variants in a given genome. Here we report three infections in southern France with a Delta 21J_AY.4-Omicron 21K/BA.1 "Deltamicron" recombinant. The hybrid genome harbors signature mutations of the two lineages, supported by a mean sequencing depth of 1163-1421 reads and a mean nucleotide diversity of 0.1%-0.6%. It is composed of the near full-length spike gene (from codons 156-179) of an Omicron 21K/BA.1 variant in a Delta 21J/AY.4 lineage backbone. Importantly, we cultured an isolate of this recombinant and sequenced its genome. It was observed by scanning electron microscopy. As it is misidentified with current variant screening quantitative polymerase chain reaction (qPCR), we designed and implemented for routine diagnosis a specific duplex qPCR. Finally, structural analysis of the recombinant spike suggested its hybrid content could optimize viral binding to the host cell membrane. These findings prompt further studies of the virological, epidemiological, and clinical features of this recombinant.

Characteristics of the first 1119 SARS-CoV-2 Omicron variant cases, in Marseille, France, November-December 2021.

One thousand one hundred and nineteen cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant cases have been diagnosed at the Institut Hospitalo-Universitaire Méditerranée Infection, Marseille, France, between November 28, 2021, and December 31, 2021. Among the 825 patients with known vaccination status, 383 (46.4%) were vaccinated, of whom 91.9% had received at least two doses of the vaccine. Interestingly, 26.3% of cases developed SARS-CoV-2 infection within 21 days following the last dose of vaccine suggesting possible early production of anti-SARS-CoV-2 facilitating antibodies. Twenty-one patients have been hospitalized, one patient required intensive care, and another patient who received a vaccine booster dose died.

SARS-CoV-2 reinfection and COVID-19 severity.

SARS-CoV-2 reinfection rate is low. The relative severity of the first and second episodes of infection remains poorly studied. In this study, we aimed at assessing the frequency of SARS-CoV-2 reinfections and comparing the severity of the first and second episodes of infection. We retrospectively included patients with SARS-CoV-2 positive RT-PCR at least 90 days after clinical recovery from a COVID-19 episode and with at least one negative RT-PCR after the first infection. Whole genome sequencing and variant-specific RT-PCR were performed and clinical symptoms and severity of infection were retrospectively documented from medical files. A total of 209 COVID-19 reinfected patients were identified, accounting for 0.4% of positive cases diagnosed from 19 March 2020 to 24 August 2021. Serology was performed in 64 patients, of whom 39 (60.1%) had antibodies against SARS-CoV-2 when sampled at the early stage of their second infection. Only seven patients (3.4%) were infected twice with the same variant. We observed no differences in clinical presentation, hospitalization rate, and transfer to ICU when comparing the two episodes of infections. Our results suggest that the severity of the second episode of COVID-19 is in the same range as that of the first infection, including patients with antibodies.

Introduction of the SARS-CoV-2 Beta variant from Comoros into the Marseille geographical area.

BACKGROUND: We describe the epidemiology of the first cases diagnosed in our institute of infections with the SARS-CoV-2 Beta variant and how this variant was imported to Marseille. METHODS: The Beta variant was identified based on analyses of sequences of viral genomes or of a spike gene fragment obtained by next-generation sequencing using Illumina technology, or by a real-time reverse-transcription-PCR (qPCR) specific of the Beta variant. RESULTS: The first patient diagnosed as infected with the SARS-CoV-2 Beta variant was sampled on January 15, 2021. Twenty-nine patients were diagnosed in January 2021 (two weeks). Fifteen (52%) patients were of Comorian nationality. Eight (28%) had travelled abroad, including six who had returned from Comoros. Phylogeny based on SARS-CoV-2 genomes from 11 of these patients and their best BLAST hits from the GISAID database showed that seven patients, including the four returning from Comoros, were clustered with 27 other genomes from GISAID that included the six first Beta variant genomes described in Comoros in January 2021. CONCLUSIONS: Our analyses highlight that, as for the case of other SARS-CoV-2 variants that have been diagnosed in Marseille, the Beta variant was imported to Marseille through travel from abroad. It had limited spread in our geographical area.

SARS-CoV-2 Vaccination and Protection Against Clinical Disease: A Retrospective Study, Bouches-du-Rhône District, Southern France, 2021.

From January 18th to August 13th, 2021, 13,804 unvaccinated and 1,156 patients who had received at least one COVID-19 vaccine dose were tested qPCR-positive for SARS-CoV-2 in our center. Among vaccinated patients, 949, 205 and 2 had received a single, two or three vaccine doses, respectively. Most patients (80.3%) had received the Pfizer-BioNTech vaccine. The SARS-CoV-2 variants infecting vaccinated patients varied over time, reflecting those circulating in the Marseille area, with a predominance of the Marseille-4/20A.EU2 variant from weeks 3 to 6, of the Alpha/20I variant from weeks 7 to 25, and of the Delta/21A variant from week 26. SARS-CoV-2 infection was significantly more likely to occur in the first 13 days post-vaccine injection in those who received a single dose (48.9%) than two doses (27.4%, p< 10-3). Among 161 patients considered as fully vaccinated, i.e., >14 days after the completion of the vaccinal scheme (one dose for Johnson and Johnson and two doses for Pfizer/BioNTech, Moderna and Sputnik vaccines), 10 (6.2%) required hospitalization and four (2.5%) died. Risks of complications increased with age in a nonlinear pattern, with a first breakpoint at 54, 33, and 53 years for death, transfer to ICU, and hospitalization, respectively. Among patients infected by the Delta/21A or Alpha/20I variants, partial or complete vaccination exhibited a protective effect with a risk divided by 3.1 for mortality in patients ≥ 55 years, by 2.8 for ICU transfer in patients ≥ 34 years, and by 1.8 for hospitalization in patients ≥ 54 years. Compared to partial vaccination, complete vaccination provided an even stronger protective effect, confirming effectiveness to prevent severe forms of COVID-19.

Limited spread of a rare spike E484K-harboring SARS-CoV-2 in Marseille, France.

We detected SARS-CoV-2 of PANGO lineage R.1 with the spike substitution E484K in three patients. Eleven other sequences in France and 8,831 worldwide were available from GISAID, 92% originating from Japan. The three genome sequences from our institute were phylogenetically closest to another from Guinea-Conakry, where one of the patients had travelled. These viruses did not exhibit any unusual features in cell culture. Spike structural predictions indicated a 1.3-time higher transmissibility index than for the globally spread B.1.1.7 variant but also an affinity loss for gangliosides that might have slowed dissemination. The spread of new SARS-CoV-2 mutants/variants is still not well understood and therefore difficult to predict, and this hinders implementation of effective preventive measures, including adapted vaccines.

Incidence and Outcome of Coinfections with SARS-CoV-2 and Rhinovirus.

BACKGROUND: We aimed to compare the clinical severity in patients who were coinfected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and rhinovirus or monoinfected with a single one of these viruses. METHODS: The study period ranged from 1 March 2020 to 28 February 2021 (one year). SARS-CoV-2 and other respiratory viruses were identified by real-time reverse-transcription-PCR as part of the routine work at Marseille University hospitals. Bacterial and fungal infections were detected by standard methods. Clinical data were retrospectively collected from medical files. This study was approved by the ethical committee of our institute. RESULTS: A total of 6034/15,157 (40%) tested patients were positive for at least one respiratory virus. Ninety-three (4.3%) SARS-CoV-2-infected patients were coinfected with another respiratory virus, with rhinovirus being the most frequent (62/93, 67%). Patients coinfected with SARS-CoV-2 and rhinovirus were significantly more likely to report a cough than those with SARS-CoV-2 monoinfection (62% vs. 31%; p = 0.0008). In addition, they were also significantly more likely to report dyspnea than patients with rhinovirus monoinfection (45% vs. 36%; p = 0.02). They were also more likely to be transferred to an intensive care unit and to die than patients with rhinovirus monoinfection (16% vs. 5% and 7% vs. 2%, respectively) but these differences were not statistically significant. CONCLUSIONS: A close surveillance and investigation of the co-incidence and interactions of SARS-CoV-2 and other respiratory viruses is needed. The possible higher risk of increased clinical severity in SARS-CoV-2-positive patients coinfected with rhinovirus warrants further large scale studies.

Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France.

Great concerns have been raised about SARS-CoV-2 variants over the past six months. At the end of 2020, an increasing incidence of spike substitutions Q677H/P was described in the USA, which involved six independent lineages. We searched for changes to this amino acid in the sequence database of SARS-CoV-2 genomes obtained at the IHU Méditerranée Infection (Marseille, France) from 3634 patients sampled between February 2020 and April 2021. In seven genomes (0.2%), we found a deletion of five amino acids at spike positions 675-679 (QTQTN) including Q677, and in 76 genomes (2.3%) we found a Q677H substitution. The 83 genomes were classified in ten different Pangolin lineages. Genomes with a spike Q677 deletion were obtained from respiratory samples collected in six cases between 28 March 2020 and 12 October 2020 and in one case on 1 February 2021. The Q677H substitution was found in genomes all obtained from respiratory samples collected from 19 January 2021 and were classified in seven different lineages. Most of these genomes (41 cases) were of UK variant. Two others were classified in the B.1.160 Pangolin lineage (Marseille-4 variant) which was first detected in July 2020 in our institute but was devoid of this substitution until 19 January 2021. Also, eight genomes were classified in the A.27/Marseille-501 lineage which was first detected in our institute in January 2021 and which either harboured or did not harbour the Q677H substitution. Thus, the spike Q677H substitution should be considered as another example of convergent evolution, as it is the case of spike substitutions L18F, E484K, L452R, and N501Y which also independently appeared in various lineages.